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Abstract. We show that, if solar 5 min. oscillations are 
excited by convection in the upper layers of the convec- 
tive envelope, it is impossible to explain the opposite line 
asymmetries observed in the velocity and intensity spec- 
tra with assumptions on the dissipations which reduce the 
problem to a second order one. The interpretation of that 
observation requires to solve the full non-adiabatic prob- 
lem which is of the fourth or sixth order. We also analyze 
the causes of line asymmetries in the frame of the general 
problem and we show that to locate the source, it is better 
to study line asymmetries not too far from line centers. 



1. Introduction 

Duvall et al. (1993) have discovered that lines of the solar 
acoustic power spectrum show asymmetries which have 
opposite signs in the velocity and the intensity signals. 
Following Abrams and Kumar (1996), a line will be said 
to have a positive (negative) asymmetry if it has more 
(less) power on the high frequency-side than on the low- 
frequency one. The velocity lines show a negative asym- 
metry while the intensity ones show a positive one though 
the amount of asymmetry varies with frequency. This find- 
ing has been recently confirmed by MDI which is one of 
the SOHO experiments. 

The line asymmetry was also predicted theoretically 
(Gabriel 1992, 1993, 1995) and discussed by Kumar 

(1994) , Lou and Fan (1995), Roxburgh and Vorontsov 

(1995) , Abrams and Kumar (1996), Rast and Bogdan 
(1997a) and Nigam et al. (1997) but in nearly all cases 
for the velocity signal only. 

It is considered that solar p-modes are non-linearly ex- 
cited by convection in the upper layers of the convective 
envelope. The theory was originally proposed by Goldre- 
ich and Keeley (1987) and further developed by Goldreich 
and Kumar (1988, 1990) and Goldreich et al. (1994) (see 
also Osaki (1990) and Musielak (1994)). Recently Rast 
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(1997b, 1997c) has proposed a slightly different mecha- 
nism in which excitations are associated with new down- 
flow plume formation. However, so far, these theories have 
been unable to predict the location and the width of the 
excitation zone which must be found from the interpreta- 
tion of observations. Attempts have been made by Kumar 
(1994) and Abrams and Kumar (1996). 

Until recently, no attempt had been made to use in- 
tensity line asymmetries. The reason for this situation is 
easy to understand. Velocity measurements are much eas- 
ily connected to theory as even the simple adiabatic theory 
makes predictions concerning the velocity. In the studies 
of line asymmetries done so far (with the exception of 
Kumar 1994), dissipations have however been taken into 
account but in a very rough way which allows, with the 
Cowling approximation, to keep a second order boundary 
value problem which is much easier to handle than the full 
fourth order non-adiabatic one. Even if we can question 
the validity of the approximations done to keep a second 
order problem when it comes to make theoretical predic- 
tions accurate enough for the interpretation of observa- 
tions, such predictions can be done and it is tempting to 
confront them with observations. However the perturba- 
tion of the luminosity does not appear in those simplified 
problems. Moreover the observations are often made in a 
frequency interval and it is not obvious whether the in- 
tensity fluctuations may be simply associated to the per- 
turbations of the total luminosity obtained by theoretical 
computations. Even if this is allowed, we have to solve a 
fourth order problem the theoretical basis of which are 
poorly understood as it requires to solve the linear inter- 
actions between pulsation and convection. It is probably 
for this reason that attempts trying to explain the oppo- 
site line asymmetries in velocity and intensity spectra with 
second order boundary value problems, have been made 
only recently after the confirmation of the discovery by 
MDI. To do that, it is necessary to introduce one more 
hypothesis which links the intensity fluctuations to either 
the Lagrangian or the Eulerian temperature (or pressure) 
perturbations. 

In this note we prove that the solution of second order 
problems leads always the same line asymmetries for the 
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two variables. It is therefore impossible to explain both the 
velocity and intensity observations with such simplified 
problems. Their interpretation requires to solve the full 
non-adiabatic equations. 

2. Equations of the problem 

To study the problem of pulsations excited by convec- 
tion, we have to solve a non-homogeneous linear problem 
(Gabriel 1993): 



dY 

dr 



AY + F 



The homogeneous system is given by the equations of non- 
adiabatic stellar stability and F is the forcing term pro- 
duced by convection. A is a matrix, function of the dis- 
tance to the center and of the oscillation frequency which 
indeed takes real values only (while the eigenvalues of the 
homogeneous problem are complex) . The solutions of this 
system have to fulfill the same boundary conditions as the 
homogeneous one. For radial oscillations and in the Cowl- 
ing approximation for non-radial ones, the system is of 
dimension four, while when the perturbation of the po- 
tential is taken into account in the non-radial case the 
problem is of the sixth order. 

The solution of that system is given by (Gabriel 1993) 



Y(r) = / G(r,r').F(r')dr' 
Jo 



(2) 



where G(r,r') is the Green function matrix of the homo- 
geneous system. 

Let 2N be the order of the system. The homogeneous 
problem has N independent solutions verifying the bound- 
ary conditions at the center Y\, Y2, . . . Yjy and N inde- 
pendent solutions verifying the boundary conditions at the 
surface Y"jv+i, Yn+2, ■ ■ ■ Y 2 n- The corresponding funda- 
mental matrix is M = (Yi,Y 2 , ■ ■ . , Yn, Y^+i, ■ ■ ■ , Y 2 n) 
and the Green function matrix is given by: 

JV 

Ga (r, /) = - £ M lk (r)M- 1 (/) r < r' 

k=i 

2N 

= Yl M lk (r)M^(r>) r > r' (3) 

k=N+l 



Equation (2) can also be written as 



r W = mr , 

2N rr 
k=2N 

= C k (r)Y k {r) 



r \M k (r')\ dr , 



where M k is obtained by replacing the k th column of the 
fundamental matrix by F . 

If the oscillations are excited by convection, the source 
term is different from zero below the observation level only 
and the coefficients of the solutions regular at the center 
are zero. 

If the first component of Y is u;5r and the fourth one 
is 5L, the power spectra for velocity and intensity fluctu- 
ations are given by 



2JV 



(1) 



2N 



£ C {r)Ct{r)Y l3 {r)Y: k {r) 

j=N+l k=N+l 



fc=l 



with i equal to 1 and 4 respectively. Because of the sum- 
mation, it is possible that the two spectra show differ- 
ent line asymmetries. If the Cowling approximation is 
used for non-radial oscillations, it is, for instance, pos- 
sible to choose Y3 such that Sr(R) = and Y4 such 
that SL(R) = 0. Then the line asymmetries for the ve- 
locity spectrum will be given mainly by the behaviour of 
C4, while those of the intensity spectrum will be related 
mainly to the variation of C3 with frequency. 

If, as discussed above, the system is reduced to the 
second order, then eq. (4) reduces to Y = C 2 {r)Y \. There 
is no summation and if one spectrum is associated to each 
component of Y, the two spectra show indeed the same 
asymmetries. Notice that this result is independent of the 
explicit form of the Aij and therefore of the hypothesis 
made to take dissipations into account, provided that a 
second order problem is obtained. 

Since it is necessary to solve the full non-adiabatic 
problem to explain the different line asymmetries of the 
two spectra, we can also wonder whether the source depth 
obtained from the velocity data only and using a second 
order problem can be trusted. 

After discussing the problem of velocity and intensity 
line asymmetries with a second order problem and fail- 
ing to explain the opposite asymmetries indeed, Rast and 
Bogdan (1997a) have suggested that the problem can be 
solved if the two spectra have different noice levels. This 
might be the case but their remark mostly challenges ob- 
servers who ought to remove properly the noice in their 
data analysis. If after this has been done and that infor- 
mations concerning the excitation source have been clearly 
obtained, the two spectra still show opposite line asymme- 
tries then the problem will have to be studied with the full 
fourth (in the radial case and in the nonradial one with 
the Cowling approximation) or sixth order system of non- 
adiabatic stability equations. 



(4) 3. The causes of line asymmetries 

We will analyze the causes of line asymmetries in the 
Cowling approximation since it is very good for 5 min. 
p-modes. Then the system (1) is of the fourth order and 
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it is interesting to write eq. (4) as 

-«[^ 

with V given by: 

{(T12-F2 — T13F3 + T14F4), (—T21F1 + T23F3 — T24-F4), 

(T31F1 — T32-F2 + T34F4), (— T41F1 + T42-F2 — T43F3)} 

= Tji is the determinant of the two by two matrix 
obtained from (Y"i,Y"2) after suppressing lines i and j. 
(we have kept the four components of F though F3 = if 
the third line of eq. (1) is the transfer equation.) 

We will first assume that the excitation force is well 
localized and can be represented by a delta function and 
afterwards we will generalize the discussion to an extended 
source. 

First, we discuss the denominator |M(r')|. It cancels 
for each eigenvalues. Therefore, in the vicinity of one of 
them a — a R + iaj, \M(r')\ = f(r', <t)(uj — o R — ioi) with 
f(r',a) ^ and for real u, \M(r')\ 2 has minima close to 
the real part of the eigenvalues. Therefore it is nearly a 
periodic function (with a "period" equal to the frequency 
separation between two successive line centers) but not 
exactly, for two reasons: 

1 . the real parts of the eigenvalues are not exactly equidis- 
tant. 

2. the extrema of |M(r')| 2 show a variation with fre- 
quency, especially close and above the cut-off fre- 
quency. This behaviour is already seen in simple ide- 
alized problems (Gabriel 1992). 

Therefore, close to the real part of an eigenvalue, |M(r')| 2 
is symmetric but over a wider frequency range, it produces 
some skewness of the line profiles. However the asymme- 
tries introduced by this term will be the same for the ve- 
locity and for the intensity spectra. 

Let us now consider the numerator. 
The two solutions regular at the surface vary slowly with 
frequency and they have generally been considered as con- 
stant in previous discussions. This is a good approxima- 
tion for narrow lines. However, if line profiles are con- 
sidered over a frequency range of the order of the eigen- 
value separation, these two solutions will also introduce 
some skewness different for each spectrum. We now dis- 
cuss the consequences of the variation of V assuming that 
Y 3 and Y4 are constant. Since V varies with the source 
type, the asymmetries will also be a function of the source 
type. Different spectra will have the same asymmetries 
only if V does not change its direction when uu varies, i.e. 
if V(r',uj) — f(uj)V 1 (r / ). This is very unlikely to occur 
as it requires very special properties for Y\, Y2 and F 
and the velocity and intensity spectra will show different 
line asymmetries though not necessarily of opposite sign. 



Because the radial and horizontal displacements of the 
eigenfunctions have one more node when the order is in- 
creased by one, the T^- show also an oscillatory behabiour 
with a quasi-period in frequency a little larger than twice 
the frequency separation between successive eigenvalues. 
However as these terms have to be squared to get the spec- 
trum, they finally lead to a function with a quasi-period 
in frequency a little larger than the frequency separation 
between successive eigenvalues. Again, because the real 
part of the eigenvalues are not exactly equidistant and 
the extremal values of the vary with frequency, the 
numerator will produce some asymmetry in the lines if a 
large enough frequency domain is considered. But more 
important, the two quasi-periodic functions appearing in 
the numerator and denominator have different "periods" 
and they will be out of phase except, may be, for one or 
two lines. (The number of symmetric lines increases with 
the source depth and if it is close enough to the surface, 
there may be no such line.) This will be the main source 
of asymmetry close to the center of the lines and the only 
one connected to the source position. Therefore to locate 
the source, it is better to study the asymmetries in a not 
too wide frequency domain around the line centers. Also, 
if a symmetric line is observed, the asymmetry presently 
under discussion must have opposite signs for lines on each 
side of the symmetric one. 

In the case of an extended source, we must notice that 
the quasi-period of V changes with the position in the 
Sun; it increases as the point under consideration moves 
inward while that of the denominator does not change and 
therefore the phase lag between the numerator and the de- 
nominator changes too. Summing up the contributions of 
several layers will indeed smooth out the oscillations pro- 
duced by the numerator and therefore reduce the amount 
of line asymmetries. 

4. Conclusions 

In this note, we have proved that the opposite line asym- 
metries observed in velocity and intensity power spectra 
cannot be explained by theories which approximate the 
dissipations in such a way to reduce the problem to a sec- 
ond order one. The solution of that problem requires the 
resolution of the full non-adiabatic problem. We have also 
analyzed the causes of line asymmetries in the frame of 
the general problem. The solution of that problem should 
be obtained including all the interactions between pulsa- 
tion and convection (i.e. not just the perturbation of the 
convective luminosity) because they are all important in 
the upper layers of the convective enveloppe. There are 
several causes of asymmetry but only one of them can be 
related to the properties of the source. To pick up that 
one it is better to study line asymmetries not too far from 
line centers. 

References 

Abrams, D. & Kumar, P, 1996, ApJ, 472, 882 



4 



M. Gabriel: On Velocity and Intensity Line Asymmetries 



Duvall, T.L.,Jr., Jefferies, S.M., Harvey, J.W., Osaki, Y. & 

Pomeramtz, M.A., 1993, ApJ 410, 829 
Gabriel, M., 1992, A&A 265, 771 
Gabriel, M., 1993, A&A 274, 935 
Gabriel, M., 1995, A&A 299, 245 
Goldreich, P. & Keeley, D.A., 1977, ApJ 212, 243 
Goldreich, P. & Kumar, P., 1988, ApJ 326, 462 
Goldreich, P. & Kumar, P., 1990, ApJ 363, 694 
Goldreich, P., Murray, N. & Kumar, P., 1994, ApJ 424, 466 
Kumar, P., 1994, Ap.J., 428,827 
Lou, Y.-Q. & Fan, Z.-H., 1995, A&A, 299, 703 
Musielak, Z.E., Rosner, R., Stein, R.F. & Ulmschneider, P, 

1994, ApJ 423, 474 
Nigam, R., Kosovichev, A.G. & Scherrer, P.H., 1997 to appear 

in Sounding Solar and Stellar Interior 
Osaki, Y., 1990, Progress of Seismology of the Sun and Stars, 

ed., Y. Osaki & Shibahashi, H., Springer- Verlag, pg. 75 
Rast, M.P. & Bogdan, T.J., 1997a, preprint 
Rast, M.P., 1997b, SCORE'96: Solar Convection and 

Oscillations and their Relationship, ed. F.P. Pijpers, 

J. Christensen-Dalsgaard & OS. Rosenthal, in press 
Rast, M.P., 1997c, preprint 

Roxburgh, I, W. & Vorontsov, S.V., 1995, MNRAS, 272, 850 



This article was processed by the author using Springer- Verlag 
M? E X A&A style file L-AA version 3. 



